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ABSTRACT 

The five enigmatic "Cocoon stars" after which the Quintuplet cluster was christened have 
puzzled astronomers since their discovery. Their extraordinary cool, featureless thermal spectra 
have been attributed to various stellar types from young to highly evolved, while their absolute 
luminosities places them among the supergiants. We present diffraction-limited images from the 
Keck 1 telescope which resolves this debate with the discovery of rotating spiral plumes charac- 
teristic of colliding-wind binary "pinwhcel" nebulae. Such elegant spiral structures, found around 
high- luminosity Wolf-Rayet stars, have recently been implicated in the behavior of supernovae 
lightcurves in the radio and optical. 



The five enigmatic "Cocoon stars" after which 
the Quintuplet cluster was christened have puz- 
zled astronomers since their discovery (8). Hun- 
dreds of stars have now been identified within the 
cluster (1, 2), placing it among the most massive 
in our Galaxy, yet the nature of the five extremely 
red stars at the heart of the Quintuplet has re- 
mained elusive. Their extraordinary cool, feature- 
less thermal spectra (~780-1315K (1)) have been 
attributed to various stellar types from young to 
highly evolved, while their absolute luminosities 
places them among the supergiants (10 4_5 L Q ). 
We present diffraction-limited images from the 
Keck 1 telescope which resolves this debate with 
the discovery of rotating spiral plumes character- 
istic of colliding-wind binary "pinwheel" nebulae. 
These have previously been reported in dust shells 
around luminous hot Wolf-Rayet stars (11, 6). 

Using high-resolution speckle techniques in the 
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near-infrared (13), all five Cocoon stars were (at 
least partially) resolved, and images recovered of 
the two presenting the largest apparent size, Q 2 
and Q3, (Fig. 1). Outflow plumes depicted follow 
the form of an Archimedian spiral, thereby estab- 
lishing the presence of circumstellar dust formed 
in a colliding-wind binary system. These rare Pin- 
wheel nebulae result when dust condensation in 
the stellar wind is mediated by the presence of a 
companion star, as established in the prototype 
systems WR 104 (11) and WR 98a (6). Dust nu- 
cleation is enabled by the wind compression associ- 
ated with the bowshock between the stellar winds. 
Newly-formed hot dust streaming into the wake 
behind the companion is wrapped into a spiral by 
the orbital motion as it is embedded within the 
expanding Wolf-Rayet wind. 

With high-resolution images available from two 
epochs separated by 357 days, the dust plume of 
star Q 3 was fitted to an Archimedian spiral model 
with winding angle 110±10 mas/turn and with an 
inclination of the rotation axis to the line-of-sight 
of 26° (not well constrained; adequate fits possible 
in the range 0-36°). The proper motion of struc- 
tures between the two epochs indicates a rotation 
period of the spiral, and hence the colliding-wind 
binary, of 850±100days. Assuming an 8000 pc 
distance to the Quintuplet (9), these measure- 
ments constrain the Wolf-Rayet wind velocity to 



be Voo = 1800±300 km/sec, which is typical for a 
late- type Carbon-rich (WC-spectrum) star (4). 

Simple models fall short of reproducing all 
structures (Fig. 1), particularly near the bright 
core where multiple knots and streams can exist. 
This has also been noted in WR 98a (6) where 
extra complexity was attributed to optical depth 
and line-of-sight effects from a 3-D structure (3). 
Q 2 appears to have similar parameters to Q 3 al- 
though a second image epoch is required to mea- 
sure the rotation rate. Partially-resolved objects 
Ql, Q4 and Q9 (see 13 and Table SI for obser- 
vational results) exhibit similar colors and surface 
brightness to Q2 & Q3, and we therefore sug- 
gest that they are also pinwheels but with tighter 
winding angles (therefore shorter periods), or less 
favorable inclinations. Furthermore, the proto- 
type pinwhcel WR 104, with a period of 243 days, 
would give an apparent size in close accord with 
measured sizes of Ql, Q4 and Q9 if it were re- 
moved to the distance of the Quintuplet. 

Given the extreme visible extinction (A v =29±5 
{2)), small separation of the central binary stars 
(likely ^0.6 mas or 5AU), and presence of high- 
luminosity circumstcllar dust shells, it would be 
extremely difficult to detect or study these sys- 
tems with other techniques. The most luminous 
stars in our Galaxy are often surrounded by dusty 
shells, and the implication that most, if not all 
of these harbor massive binaries (not single stars) 
has important ramifications for the high-mass tail 
of the stellar Initial Mass Function (IMF). Bina- 
rity is also a key element to studies of Type Ib/c 
and Type lib supernovae. There are recent indi- 
cations that explosion lightcurves can be modified 
by the imprint of circumstcllar matter, carrying an 
encoding of the mass-loss history of the supernova 
precursor star system (7, 10). 
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Fig. 1. — False-color images of Q 2 at 3.08 /mi 1999 July (left panel) and Q3 at 2.21 /mi from 1998 Aug 
(upper right panel) and 3.08 /mi from 1999 July (lower right panel). Overplotted on the Q3 images is a 
rotating Archimedian spiral model fitted to the dominant tail of the outflow plume at the two separate epochs 
(dashed line). Identification of Quintuplet objects including Q2, Q3 is discussed further in the supporting 
online material and Fig. SI. 
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Supporting On-line Material 

1. Observations and Discussion 

Observations of the Quintuplet cluster are 
listed in Table SI and were made with the near- 
infrared NIRC camera on the Keck 1 telescope. 
The nomenclature of Moncti et al.(5) is adopted 
here, identifying the dusty red Quintuplet proper 
members, or cocoon stars, as Ql, Q2, Q3, Q4 
and Q 9. Locations of each of these stars within 
the cluster as imaged by the Hubble Space Tele- 
scope is given in Fig. SI, which also depicts our 
diffraction-limited Keck images of Q2, Q3 over- 
layed with graphical indication of the relative spa- 
tial scale. 

The high angular resolution imaging utilized 
rapid-exposure speckle interferometry techniques 
from a number of long-standing experiments at the 
Keck 1 (8). Typical datasets entailed from one to a 
few hundred short-exposures (~0.14sec) with the 
NIRC camera operating in a high-magnification 
(0.0206 arcsec/pixel) mode. Interleaved observa- 
tions of galactic-center object IRS 7 were used 
to calibrate the telescope-atmosphere point-spread 
function (the unresolved nature of this object 
was itself checked against stars HD 159255 and 
HD 163042). 

All five Quintuplet cocoon stars were found to 
be spatially resolved to some degree. The full- 
width at half-maximum (FWHM) of best-fitting 
circular Gaussian functions, used here to give 
an estimate of apparent size, are given in Ta- 
ble SI. Images of stars Q2 and Q3, recovered 
with Fourier techniques (8) and presented earlier in 
Fig. 1, show that complex and asymmetric struc- 
tures exist within the dust shells surrounding these 
stars. The complexity of these visibility functions 
resulted in poor fits for the Gaussian model, and 
we therefore limited the fit range to low spatial 
frequencies/short baselines where visibility excur- 
sions are not yet strongly manifest (<2xl0 6 and 
lxl0 6 rad _1 for 2.21 and 3.08 /im respectively: 
equivalent to ^4.5 and 3.1m). 

By implication, if similar physics is assumed to 
pertain to the remaining three of the five Quintu- 
plet cocoon stars, then the utility of fitting Gaus- 
sian shells may appear limited. Despite this, for 
objects which are only partially resolved, the fit- 
ting of such simple functional forms to interfer- 



ometer data gives a good estimate of the over- 
all size, and permits quantitative comparison with 
models. Furthermore, this allows sizes for mor- 
phologically complex objects Q 2 and Q 3 to be 
compared against measurements of the partially- 
resolved Q 1, Q 4 and Q 9, and against other dusty 
Wolf-Rayets observed with high resolution tech- 
niques. Note also that FWHM in the range 10- 
15 mas are only marginally resolved with this ex- 
periment; relative errors in these cases are corre- 
spondingly high. 

Apparent sizes have been combined with flux 
measurements to yield estimates of surface bright- 
ness in our two filter bandpasses for the target 
stars. Correction for the A v =29±5 visible extinc- 
tion was made using the optical dust constants 
of Mathis(4). Surface brightnesses were then de- 
rived following Monnicr ct al.(6), which also gives 
a comparison population of Galactic dusty Wolf- 
Rayets with similarly measured surface bright- 
nesses. For both filter bandpasses, the Quintuplet 
cocoon stars showed surface brightnesses within 
the range spanned by the WR population studied 
in Monnier et al.(6). In particular, for Q4 and Q9 
the color temperature between 2.21 and 3.08 ^m 
appeared to be in reasonable accord with similar 
measurements from the Galactic population. This 
finding was found to be generally robust against 
variations in the extinction correction over the ex- 
pected range. 

However, in some regards the Quintuplet WRs, 
and in particular Q2 and Q3, did appear to be 
distinct from the Galactic population. The ob- 
served increase in size between 2.21 and 3.08 /im 
which approaches a factor of ~2 in the Quintuplet, 
is found to be a more modest ~1.4 elsewhere(9; 6). 
Such changes of size (and surface brightness) with 
wavelength reflect the fact that dust with a range 
of temperatures contributes to the near-IR emis- 
sion. The dramatic enlargement between 2.21 and 
3.08 /im argues for a flatter thermal profile in the 
Quintuplet dust shells: plausibly due to external 
heating from stars in the dense central region of 
the cluster. The comparison population of Galac- 
tic WRs were in far less crowded regions where the 
outer dust shell likely receives little or no energy 
except that originating with the central WR star. 

Perhaps the simplest way to test the hypoth- 
esis that the remaining 3 cocoon stars are also 
pinwheels is by photometric monitoring. Variabil- 
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ity of these sources has already been noted(2), 
although not with sufficient coverage to reveal 
a cyclic change consummate with any rotational 
period(3). The ^2.5 yr periods inferred from the 
imaging extend the confirmed operation of contin- 
uous dust formation by the pinwheel mechanism 
to significantly longer periods/larger binary sepa- 
rations than previously known, with implications 
for models of these processes. 

This finding of late-type WC binaries in the 
Quintuplet means that WC (Carbon rich) out- 
number WN (Nitrogen rich) Wolf-Rayets by 11:6, 
and furthermore all WC stars are dusty. This 
makes an interesting contrast to the massive young 
WR-rich cluster Westerlund 1 where the WC:WN 
ratio is reversed to 7:12, with none of the WC's 
exhibiting dust. Clearly, the close binaries at the 
heart of the Pinwheel systems can be responsi- 
ble for significant modification of the stellar evolu- 
tionary path, in which mass-transfer or envelope- 
stripping events might precipitate the WC phase. 
This entanglement of binarity, mass-loss history, 
and evolutionary path has the potential to skew 
population distributions, although it is unclear ex- 
actly what conditions resulted in the abundance of 
Pinwhcels in this cluster. 
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Table SI 
Observing log and apparent sizes 



Source(5) 


Alt. Name (7) 


Date 


2.21 /zm FWHM (mas) 


3.08 A*m FWHM (mas) 


Ql 


GCS 3-4 


1998 Aug 06 


18±3 




Q2 


GCS 3-2 


1998 Aug 06 


35±2 








1999 May 04 


38±3 








1999 Jul 29 


37±3 


77±2 


Q3 


GCS 4 


1998 Aug 06 


41±2 








1999 Jul 29 


40±2 


75±2 


Q4 


GCS 3-1 


1998 Aug 06 


13±4 








2002 Jul 23 


15±3 


20±3 


Q9 


GCS 3-3 


1998 Aug 06 


13±4 








2002 Jul 23 


< 11 


21±3 



Note. — Log of observations of the five cocoon stars. Observations in the K band at 2.21 /im 
were made at all epochs, but longer wavelength data at 3.08 /im were only secured on the four 
occasions listed. The FWHM of a circular Gaussian profile fit to the visibility data is also given, 
together with the estimated uncertainty, as a measure of the overall apparent size. 
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Fig. SI- The background star-field is from multi- 
wavelength Hubble Space Telescope (NICMOS) 
near-infrared imaging. Further details and discus- 
sion of this image can be found in Figer et al. 
(1). The five dusty red cocoon stars are labelled 
according to the nomenclature of Moneti et al. 
(5). Inset images of Q2 and Q3 recovered with 
our Keck imaging experiments (see also Fig. 1) 
are overlayed, with graphical indication showing 
the relative scaling between the Hubble and Keck 
imaging. 



